Zero energy Majorana modes in superconducting wires. 
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I consider a model of two thin superconducting wires interacting through Josephson tunneling 
and Coulomb interactions. It is shown that there is a parameter range where such a model may 
have two ground states with different types of quasi-long-range order: one is superconducting and 
the other is pair density wave. It is further shown that boundaries between these different phases 
support zero energy Majorana modes and that these modes are robust against disorder. 



Majorana fermions have attracted a great deal of at- 
tention lately. A part of this interest is related to the fact 
that zero energy Majorana modes (ZEMMs) are nonlo- 
cal and robust with respect to all sorts of perturbations. 
The latter property is related to the fact that ZEMMs 
emerge as bound states of Majorana fermions and topo- 
logical objects such as vortices [H , [H , [H or kinks [4]. The 
nonlocality is related to the fact that to define an occu- 
pation number operator one needs at least two ZEMMs. 
More formally, the dimensionality of the Hilbert space 
spanned by N degenerate Majorana modes is 2^ N / 2 \ In 
particular, to make a two-level system one needs at least 
two modes located at different positions. Thus spatially 
separated ZEMMs correspond to the situation described 
in the Einstein-Podolsky-Rosen paradox. 

The Majorana fermion is not as exotic as it might seem, 
being a linear superposition of an electron and a hole, 
and hence a particular case of a Bogolyubov quasipar- 
ticle. Indeed, the most straightforward way of realizing 
Majorana fermions is to create them literally from elec- 
trons and holes and they have been suggested to emerge 
in p-wave superconductors [I[,[4[, topological insulators 
(see, for instance 0,0) or in semiconducting wires with 
strong spin-orbit coupling in contact with s-wave super- 
conductors 7\. There are other possibilities, however. 
Majorana fermions may emerge as collective excitations 
as, for instance, in the one-dimensional transverse field 
Ising model or in the Kitaev model 0. In these cases 
Majorana fermion operators are constructed from spin 
S=l/2 ones via a Jordan- Wigner transformation; the re- 
sulting objects are highly nonlocal in spins. 

In the current paper I consider Majorana fermions in 
one-dimensional superconductors. However, instead of 
being Bogolyubov quasiparticles, they emerge as collec- 
tive excitations; more precisely, as nonlocal objects made 
of superconducting pairs. Such a route may appear to be 
unnecessarily complicated, but, in fact, represents several 
advantages. One of these is that to get the desired re- 
sult one does not need any exotic superconductors, like p- 
wave ones, but can manufacture Majorana fermions from 
conventional s-wave or high temperature d-wave super- 
conductors. Here I am encouraged by the experimen- 
tal observation of the magnetic field driven insulator- 
superconductor crossover in superconducting wires [9( 
and by its theoretical explanation given in [10| . 

Below I consider a model of two coupled superconduct- 



ing wires. As is known from previous work (Io| . such a 
system can be in two different phases depending on the 
balance between the inter-wire Coulomb repulsion and 
Josephson tunneling energy. One phase supports super- 
conducting quasi-long-range order (SC) and the other one 
supports pair density wave (PDW) quasi-long-range or- 
der. In the presence of disorder the PDW is pinned and 
this phase is an insulator. I demonstrate that boundaries 
between the SC and PDW phases support ZEMMs and 
this fact is not affected by disorder. I also show that 
in the state where domain walls between SC and PDW 
phases have a finite density a narrow Majorana band is 
formed with a density of states which has a sharp singu- 
larity of zero energy. 

Consider two coupled superconducting wires separated 
by a thin insulating layer of width w as has been consid- 
ered in [lj,[ri|,[12|. This construction can be also re- 
alized as a double chain Josephson junction array or, in 
the limiting case, as a system of two chains of attractive 
Hubbard models. Since all spin degrees of freedom are 
gapped, at low temperatures the dynamics of individual 
wires are described by the phase fields $i (i = 1,2). The 
wires are coupled by the Josephson pair tunneling and 
Coulomb interaction. 



H = H x + H 2 + V, 
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V = I dx - gj cos($i - $ 2 - Qx) + Unxn 2 , (3) 



The operators rij denote density fluctuations of the 
Cooper pairs and can be represented as |13| : 
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where ©i is the field dual to [d x Q(x), $(y)} = 

—mS(x — y). As it was argued in jlOj, Q is determined by 
the effective magnetic field and is given by the deviation 
of the vortex density from the nearest integer value N: 
Q = 2ewB/c - N. Substituting (0} into © and defining 
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the fields 



*1.2 = 



^(^ /2 <f> + ±A^ 2 $_), 



® X z = ^L(k- 1 / 2 q + ±k-_ 1/2 &- 

[d x e a {x), $ b (y)] = -iS ab 5(x - y), 



(5) 



I obtain the Hamiltonian H = H + + ii_ , where H + de- 
scribes the symmetric mode (+): 



H+ = yJ dx[(d x <S>+) 2 + (d x O + ) 2 
and iJ_ contains only the asymmetric fields: 
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cos L/4nKZ$- - Qx^j - Vj cos ^Air/K^O^ |, 



where V c = n^U, Vj — n 2 gj and 
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Depending on which of the cosines in takes over, 
the ground state of this model describes either quasi long 
range superconducting order or pair density wave. The 
latter state has a singularity in the density-density cor- 
relation function at the finite wave vector 27m . When 
both cosines are relevant (that is at 1/2 < K- < 2) these 
states are separated by a quantum critical point (QCP), 
the location of which is approximately determined by the 
relation (Vj/A) K - ~ {V c /K) 1/K -, where A is the ultra- 
violet cut-off. Fulfillment of the above condition on A_ 
is essential for the subsequent arguments. 

The vicinity of the QCP can be studied analytically 
when ss 1 (this will be my assumption throughout 
the rest of the paper). In that case it is convenient to 
refermionize ([7]) with the result 

H- = 



dx{ — (-p R d x p R + p L d x p L - rj R d x ?] R + ij L d x ri L ) H 
47ro_(A_ - ^prPlVrVl + 2im + p R p L + 2im-rj R r] L }, 



m± = Vj± V c . 
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Here pl,r and t)l,r are left- and right-moving compo- 
nents of Majorana (real) fermions. This model is equiv- 
alent to the continuum limit of two quantum Ising (QI) 
models coupled by the energy density operators. It is use- 
ful to express important operators of the original problem 
((T|) in terms of operators of the QI model. To this end I 
will use the standard correspondence between the bosonic 
exponents and the order er 1|2 and disorder pi 2 operators 
of the QI models about which all the relevant informa- 
tion can be found, for instance, in (l4j . For instance, for 



the first chain the superconducting order parameter and 
the most relevant part of the oscillatory density Q are 
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(10) 



CT1CT2 + 1M1/X2 



J210! _ 



V1P2 



(11) 
ipi<j 2 . 



Since (a) 7^ when the corresponding fermionic mass is 
positive and zero when it is negative (for (p) it is the 
other way around) , we see that in the phase where m + = 
Vj + V c > 0, m_ = Vj — V c > it is the superconducting 
order parameter which acquires a finite amplitude and 
in the phase m + > 0, m_ < this happens for the 
oscillatory part of the density. 

Neglecting the four-fermion interaction, I obtain from 
([9]) that the amplitude of the superconducting order pa- 
rameter diminishes when Vj approaches V c from above 
and vanishes at Vj < V c . Elementary scaling considera- 
tions yield the following estimate for the amplitude: 



{a 1( j 2 )^{V]-V 2 f 



(12) 



At V c > Vj one gets the same expression for the ampli- 
tude of the pair density wave. 

The transition between two phases can be driven by a 
magnetic field. A magnetic field adds to © the term 



2ivQ(r) L p L - iuPr) 



(13) 



(I drop the subscript "minus" on the velocity) . Then the 
spectrum becomes 



w?,2 = (vk) 2 + (vQ) 2 + (V 2 + V 2 ) ± 
[v 2 kQf + v 2 Q 2 Vf + V 2 Vj] 1/2 . 

The QCP where id- (k = 0) = is achieved for 

(vQ) 2 = m+m_ = (V 2 - V 2 ). 
The magnetic susceptibility 

,d 2 F, 



(14) 



(15) 



X (B = 0) = -(2ew/c) 
1 / / ^2 m+m- 



Y ln(m+/|m_|). 



(16) 



changes its sign at the critical point m_ = 0, being dia- 
magnetic in the SC state. 

At small wave vectors we have 



{vkf 



{vQf 



Vj^V 2 



(17) 

+ {Vj - VWTQ 2 ) 2 . 



Therefore the magnetic field shifts the critical point. Ac- 
cording to 10 . such a transition (at finite temperatures 
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it becomes a crossover) from SC to PDW state, insti- 
gated by magnetic field, has already been observed @- 
However, the shift of the transition point is not the only 
effect of the field. At 



(vQ) 2 > V 2 /2 + Vjyfvj/4 + V?, 



(18) 



the spectrum changes its form and instead of having one 
minimum at k — acquires two minima at 
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(vQ) 4 - [vQfV] 
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The spectral gaps are A = V c \J\ — (Vj/vQ) 2 . This leads 
to appearance of incommensurate oscillations in correla- 
tion functions. 

Now imagine that Vj, V c are coordinate dependent and 
at some point one of the Majorana masses changes sign 
(since the Josephson tunneling gj strongly depends on 
the interchain distance, this can be easily achieved by 
changing the distance between the wires). Then there is 
a zero energy Majorana mode bound to such a kink; its 
wave function is given by 



4( X ) 
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where 7 is the operator of ZEMM {j 2 = 1), W is the 
normalization factor and the dots stand for modes with 
non-zero energy. The sign in (1201) depends on whether 
m s (x) is a kink or an anti-kink. ZEMM is robust with 
respect to changes in K _ , so far as both cosine operators 
in (J7J remain relevant (1/2 < K- < 2). In a system with 
periodic boundary conditions there is always an anti-kink 
carrying another ZEMM "f(x C 2). Together two gamma- 
matrices create a two-dimensional irreducible represen- 
tation, where creation and annihilation operators are de- 
fined as j(x c {) ± i 7 (a; C 2). 

ZEMMs can be used to assist in a coherent hopping 
"teleportation" ) of massive Majorana fermions between 
different pairs of wires, as takes place in the Kitaev 
model. Consider two such pairs where the phase bound- 
aries between the SC and the PDW phases are brought 
sufficiently close to each other, as on Fig.l. ZEMMs are 
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FIG. 1: This picture illustrates a coherent tunneling ("tele- 
portation") scheme for massive Majorana fermions described 
in the text. Solid lines represent superconducting wires. Ar- 
eas with different filling are PDW and SC phases. ZEMMs 
are located on the phase boundaries and are represented by 
the ellipses. 



involved in the smooth part of the operator 

1 --d x e- 



[Til — n2)smooth — ' 

V nK - 

itp {x - x c )~f(x c )[p R (x) ± Pl(x)] 



(21) 



where ipo is the ZEMM wave function (|20l) and the sign 
depends on whether the SC phase is on the left or on the 
right. The contribution of (|2"Tj) to the repulsive interac- 
tion of pairs is 



(19) V = V(r 12 ) 7 (r 1 ) 7 (r 2 )^(r 1 ) - p L (ri)}[p R (T 2 ) + p L (r 2 )l22) 



(the signs correspond to the configuration depicted on 
Fig. 1), where V is the matrix element of the interaction 
and ri^ are located on different pairs of wires. If the 
given pair of domain walls is sufficiently isolated from 
all others so that the given ZEMMs do not interact with 
other ones, then N — 17^1)7^2) = ±1/2 is an inte- 
gral of motion and can be replaced by a c-number. Thus 
the four-fermion term (|22[) becomes a two-fermion term 
describing a coherent hopping of massive p-fermions be- 
tween the pairs of wires. 

It is interesting to note that the expression for the total 
superconducting current 



J 



ep s 
2mc 



ey/2TrK+p s 

dx^2) = — ^ o x <& + , 

2mc 



(23) 



does not depend on the masses. So, provided the stiff- 
ness p s does not change through the system, the current 
is the same in the SC and PDW phases! This reflects the 
fact that in the absence of disorder the PDW state sup- 
ports a sliding density wave which is equivalent to phase 
coherent tunneling of superconducting pairs through the 
PDW region. However, with disorder present, the phase 
6 + is likely to get pinned so that the PDW regions be- 
come insulating. In the present context this was noted 
in 10]. It should be added however, that the Josephson 
coupling between two SC regions separated by a PDW 
one now depends on the amount of disorder in the PDW 
region and can, at least in principle, be made as strong 
as one likes by cleaning the sample. 

From the point of view of quantum information appli- 
cations such as considered in [15[ the most interesting sit- 
uation is when boundaries between SC and PDW phases 
are few and far between. Then these modes do not in- 
teract and what remains of them is the Clifford algebra 
of 71,. ..77V operators (N is the number of boundaries) 
which is equivalent to the algebra of Dirac 7-matrices. 
However, the situation when the zero modes overlap and 
form an impurity band is also interesting. This situation 
was considered in [16|,ll7|, where it was found that the 
band contains one delocalized state at zero energy and 
gives rise to a singular contribution to the specific heat: 



C(T)~ ni ln- 3 (T /T), 



(24) 



where rij is the concentration of the boundaries and To is 
a nonuniversal energy scale of the order of the impurity 
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band width. As far as I can see, ZEMMs do not yield sin- 
gular contributions to any other thermodynamic quanti- 
ties or correlation functions except the thermal conduc- 
tivity. The reason for this is that in all observables except 
the energy density ZEMMs appear in conjunction with 
operators for which correlation functions decay exponen- 
tially. 

The disordered phase described in [l6j],[I3] is some- 
what special. In the present context it is just a random 
sequence of SC and PDW phases with each the phases 
being essentially clean. In real wires there is another 
type of disorder, in the form of impurity potentials cou- 
pled to the particle density (0]). As was rightly pointed 
out in 10], this disorder becomes effective in the PDW 
phase transforming it to an insulator. I suggest, however, 
that the conventional disorder does not affect the ZEMM 
so that the results of [Hj],[l3] including formula (|24|) are 
robust. Indeed, the most relevant part of the impurity 
contribution to the Hamiltonian comes from the coupling 
of the disorder potential to the 27rno Fourier harmonics 
of the density operator: 



U (x)p(x)dx 



(25) 



da; 



U+(x)atH2 + iE/-(x)//i<7 2 



h.c 



}■ 



where £/+,?/_ are Fourier envelopes of the disorder po- 
tential with wave vectors concentrated around 2nno. In 
the SC phase (u a /Xb) = and the disorder is not effective. 
In the PDW state one of the amplitudes acquires an av- 
erage vacuum value and the most relevant contribution 
to the Hamiltonian dU) becomes 



SH = A dx\U(x)\ cos^ir/K+O+ix) - a(x)], (26) 



where A ~ (V c 2 — V]) 1 ^ and a{x) is a random phase. 
This perturbation is relevant at K + > 1/6 and leads to a 
pinning of the phase field Q+. However, the interaction 
predominantly affects the 0+ field and not the Majorana 
modes. That the interaction between the latter and + is 
weak can be ascertained from the fact that the staggered 
density amplitude vanishes at the boundary between the 
SC and PDW phases where the zero modes are located. 



Using Eqs.(83,84) from [16[ we get (x < is in the PDW 
phase): 

A(x) = (27) 
(ctijU 2 ) = (aoTO + 77j_) 1/8 exp{-/[TO_(x- x c )]}, 



#o(M) + a--i(M) 



(clq is the lattice cut-off), so that at the boundary where 
the Majorana zero mode is located the amplitude of the 
staggered density vanishes: A(x) — (ai(i2(x)) ~ \x — 

Thus I have demonstrated that in the model of cou- 
pled superconducting wires there is a parameter region 
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FIG. 2: e = A(x)/A(— oo) as a function of 2m-(x — x c ). 



where one can expect the existence of zero energy Ma- 
jorana modes. I think that in the light of experiments 
Q material realizations of such systems are entirely re- 
alistic. Local changes of the model parameters required 
to create ZEMMs can be achieved either by local mag- 
netic fields or by a change of the stripe width w. This 
may help to produce a material realization of the Kitaev 
model using conventional superconductors. 
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